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Lattice fermion actions are constructed with path integrals which are equivalent to the free one- 
flavour staggered fermion determinant. The Dirac operators used are local and have an identical 
spectrum of states to the staggered theory. Operators obeying a generalised Ginsparg- Wilson rela- 
tion are developed. 
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I. INTRODUCTION 

A complete study of QCD on the lattice requires 
the numerical simulation of dynamical fcrmions. These 
Monte Carlo calculations are extremely computationally 
costly, since the effects of quarks must be included by first 
integrating out the fermion path integral, and then de- 
scribing the resulting non-local dynamics of the fermion 
determinant. For a review of recent developments, see 
e.g. P, Q. With present techniques, the most cost- 
effective means of performing these simulations is to use 
the staggered fermion formulation of Kogut and Susskind 
[3- Recent calculations by the MILC collaboration Q| 
have demonstrated good agreement between experimen- 
tally known strong-interaction measurements and stag- 
gered fermion lattice QCD simulation. 

The formulation as it stands has a serious deficiency for 
dynamical simulations. In four dimensions, the staggered 
fermion determinant describes four flavours of fermion, 
not one. This means that while it is very simple to 
simulate four mass-degenerate fermions with the stag- 
gered method, the study of one or two flavours must 
use a fractional power of the fermion determinant. This 
raises difficult theoretical problems: what are the fermion 
fields, and what is the local action on these fields which 
reproduces this determinant? Without a path-integral 
representation of the fermion determinant, all the stan- 
dard quantum field theory construction of propagators 
(which are the two-point functions of the underlying 
quark fields) is poorly defined. If no local action exists, 
an even more severe issue arises, since there is then no 
guarantee that the continuum limit of the lattice simu- 
lation is in the same universality class as QCD and the 
link with physics is lost. 

In this paper, we describe a numerical construction of 
an operator that defines a lattice quantum field theory 
equivalent to a single, free staggered fermion. Most of 
the construction is performed in two dimensions to ease 
the computations, but some suggestive results in four di- 
mensions indicate the same construction works there too. 
Note that all the work in this paper is for the theory of 



free fermions, and the question of defining the interact- 
ing theory remains open. We do however regard this 
as a useful starting point for the more difficult problem 
of finding a path integral representation of the staggered 
fermion determinant in the presence of background gauge 
fields and the construction presented does suggest how 
to proceed further. The paper is organised as follows: 
Sec. [n] briefly describes the free staggered fermion, and 
Sec. IIIII describes the numerical construction of the local 
operator. Sec. IIVI then presents a different numerical 
construction that is seen to obey a modified Ginsparg- 
Wilson relation. In Sees. Ivl and IVTl a discussion of our 
results and conclusions is given. 



II. STAGGERED FERMIONS 

In this section, a brief overview of the staggered 
fermion formulation of Kogut and Susskind 0| is pre- 
sented, emphasising some critical properties of the re- 
sulting free quark operator. 

The staggered fermion formalism is constructed by first 
writing a naive representation of the Dirac operator on 
the lattice 
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where the Euclidean space indices (x, y) and Dirac al- 
gebra indices (i,j) have been included explicitly. This 
operator has poles not only at zero momentum, but also 
at the corners of the Brillouin zone. Counting these poles 
suggests the field coupled through this interaction matrix 
can be thought of as representing 2 d flavours of fermions 
in d dimensions. A local change of variable at every site 
of the lattice, %(x) = T(x)ip(x) with 



T(x) 



n w 



(2) 
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diagonalises the naive operator M over the Dirac algebra. 
The number of flavours is reduced by discarding all but 
one of the diagonal components of \. The operator is 
then regarded as acting on n t — 2 d /2 d / 2 — 2 d / 2 flavours 
of fermions. These flavours, all of which appear in a single 
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instance of the staggered field, are often termed "tastes" . 
The operator on the field is then 

Q x y = amS x . y + 2 X! 'fa fa) ($H-A>» - ^x-p.,y) • (3) 

where i]^(x) is the staggered phase, given by 

^(x) = (-l)£&X (4) 

A natural decomposition for the operator is to break 
the lattice into hypercubes of side-length b = 2a contain- 
ing 2 d sites. A site on the full lattice can be labelled with 
co-ordinates 

Xfl = 2Nf, + Pfl) (5) 

where are the co-ordinates of sites on the blocked 
lattice and p^ S {0,1}. The 2 d staggered variables in 



where and are the simplest representations of the 
first and second derivatives on the blocked lattice, 

A^N, N') = 4(<W,JV' ~ Sn-kN'), (10) 



□„(JV, N') = ^(6 N+fl ,N> + Sn-^n* - 25 N , N ,). (11) 

This representation makes the Dirac and taste struc- 
ture of the staggered operator more apparent. Taking 
a direct fractional power of a matrix does not change 
its structure, so it seems counterintuitive to expect the 
matrix Q 1 ^™' to be a sensible representation of the one- 
flavour Dirac operator, and this operator has been shown 
to be non-local [f|. To construct a lattice fermion with 
a more physical interpretation, begin by noting that the 
operator is 75-hermitian, namely 

=75^75, (12) 

so that 

detQ f =detQ, (13) 

and 

VdetQtQ = detQ. (14) 

The product □ = Q^Q is diagonal in the spinor index, 
and has the form 

□ = E D m- (15) 



a hypercube can be labelled by these hypercubic offset 
vectors, p 

X P (N) = X (2N + p). (6) 

Introducing a new spinor field ip ab on the blocked lattice 
sites, N 

r"(N) ^[T/ fc Xp (iV), (7) 
p 

gives the staggered fermion action in terms of these 2 d / 2 
tastes of Dirac spinors, 

S s ta g = & 4 £ $(N)Q(N, N')i>(N'), (8) 

N,N' 

with 



(9) 

I 

The lattice interaction □ thus resembles 2 d distinct copies 
of the simple discretisation of the continuum Klein- 
Gordon operator, —V 2 + m 2 on each of the 2 d lattices 
with spacing b = 2a. 



III. AN EQUIVALENT LOCAL DIRAC 
OPERATOR 

In order to define a theory with a single flavour of 
fermion, the standard method is to consider the appro- 
priate fractional powers of the fermion determinant. A 
single flavour of staggered fermion would then be repre- 
sented by det Q 1 /"*. It is important to recognise the sig- 
nificant theoretical difficulty with this prescription: the 
fractional power of the determinant can no longer be writ- 
ten directly as a path integral over Grassmann fields cou- 
pled through a local operator (Q 1 ^ nt is non-local) and 
hence all the standard quantum field theory mechanisms 
for generating correlation functions by adding sources to 
the path integral no longer follow. Locality ensures that 
interacting theories are in the same universality class of 
the continuum field theory. In order to define a sensi- 
ble lattice quantum field theory, an operator with the 
property 

det£> = detQ 1/,lt , (16) 

is required, where D defines local interactions 0- With 
this property a path integral representation can be made, 



J 



Q(N, N') = m(I ® I)8 NjN , + ^( 7m <g> I)A^N, N') + -b(j 5 ® t^D^N, A'), 
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namely 

detQ 1/nt = JV^Vip exp {-^Di/)} . (17) 

Given this form, correlation functions of the theory can 
then be constructed by adding source terms and following 
the standard construction. In this section, an operator 
D obeying Eq. 1161 is defined numerically for the free 
staggered fermion theory. 

One observation in beginning the construction of D is 
helpful: the operator will not obey the staggered-fermion 
Dirac algebra (which scatters the spin and taste compo- 
nents over the corners of the unit hypercube) since a 
counting of degrees-of-freedom suggests there are certain 
to be too many flavours. Instead, there must be sites on 
the lattice with no quark field. To construct a fermion 
field with the correct number of degrees of freedom in a 
hypercube requires there to be 2 rf / 2 components, rather 
than the 2 d components of the staggered field. To begin 
construction assume there is a single Dirac spinor at one 
site per hypercube, with no degrees of freedom on the 
other sites of the cell. 

The equivalence property of Eq. 1161 is sufficient to de- 
fine the path integral, but can be trivially satisfied for 
the free theory: any non-singular matrix can be made 
to obey this constraint after a rescaling. For the free 
case a more stringent definition of equivalence must be 
made, namely that the energy-momentum dispersion re- 
lation for fermions in the two theories be related. This 
will be satisfied if the operator itself squares to the free 
Klein-Gordon operator on the blocked lattice, i.e. if (for 
massless fermions) 

D^D = -U. (18) 

In this work, this property will be denoted "strong" 
equivalence, while the condition of Eq. EH which is trivial 
for free fermions but non-trivial in the interacting theory 
is denoted "weak" equivalence. 

The properties expected from a well defined lattice 
Dirac operator D are locality, the correct continuum 
limit for momenta below the cutoff, ir/a and invertibility 
at all non-zero momenta. Once these properties are 
satisfied the Nielsen-Ninomiya theorem |(| excludes the 
possibility of having invariance under continuous chiral 
transformations. This last issue will be dealt with in 
the next section. Here, before beginning to describe in 
detail our proposal of a Dirac operator it is helpful to 
state what locality means. An action density which has 
nearest-neighbour interactions or interactions that are 
identically zero beyond a few lattice units is certainly 
local, but no physical principle requires this extreme 



case |9J. On the lattice an action is termed local if 
its couplings have exponentially decaying tails at large 
distances. This property is ensured if D(p) is an analytic 
periodic function of the momenta p^ with period 2n/a. 

The following ansatz for a solution to the "strong" 
equivalence constraint for the blocked lattice with spac- 
ing b is made: 

D = TmP/x - q, (19) 
with p^ and q such that D obeys Eq. ^| so 

p^ -q 2 = D. (20) 

A numerical prescription for constructing an effective 
representation of the Dirac operator for massless fermions 
is used. To begin, a sequence of "ultra-local" operators of 
finite, increasing range is defined. The finite range opera- 
tor can be described with a number of coefficients weight- 
ing each distinct hopping term. The hopping terms are 
taken from A r , the set of all vectors a whose range is less 
than r. The "taxi-driver" metric is used to define the 
range of a vector a with components <2j, so 

r(a)=X>|. ( 21 ) 

i 

In two dimensions then, 
Ao = {(0,0)}, 

At = {(0,0), (1,0), (-1,0), (0,1), (0,-1)},... (22) 

A general ansatz for both p M and q, connecting fields at 
sites x and y is then 

p£ y = £ ^(a)<Wa,y, (23) 

a6A r 

and 

9x,y = ^ 9 (a)(S x+ a : y. (24) 

aeA r 

The coefficients are constrained so that the required sym- 
metries of each operator are preserved to ensure the ac- 
tion is a scalar. This implies that the coefficients u> q 
form a trivial representation of the lattice rotation group, 
while form a fundamental representation. In two di- 
mensions (where the relevant rotation groups is C± u ) the 
required irreducible representations are A\ for q and E 
for p. This in turn implies the relations 



oj q (a 1 ,a 2 ) = w g (oi, -<z 2 ) = uj q (—a\, a 2 ) = u) q (—ai, -a 2 ) = 
w 9 (a 2 , ai) = u q (a 2 ,—ai) = aj q (-a 2 , ai) = uj q (-a 2 , -ai), 



(25) 
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and 



ajp(ai,a 2 ) = Wp(ai,-a 2 
Wp(a 2 ,Oi) = -Wp(a 2 



) = -Wp(-oi,a 2 ) = -Wp(- 



l) = w p(- a 2,ai) 



a 2 , -ai). 



(26) 



One further constraint is added to improve the represen- 
tation of low-momentum states. The coefficient w g (0,0) 
is chosen such that the operator q vanishes on a zero- 
momentum plane-wave. The number of free parameters 
in the operators p and q in two and four dimensions is 
given for a few low ranges in Table [I] 
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10 
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14 


11 


5 


15 


9 


25 


17 


10 


55 


30 


189 


93 



TABLE I: The number of free parameters in the finite-range 
operators in two and four dimensional lattice actions 



A sequence of lattice Dirac operators, Di, £) 2 , ■ • ■ with 
increasing range is then considered. Each operator in 
the sequence is chosen to minimise p% , a positive-definite 
measure of the difference between the two sides of Eq. 
1181 namely 



1 



-Tr (X r 2 ), 



(27) 



Ad 2 N s 

with N s the number of sites on the blocked lattice and 
X r = D ] r D r + □. (28) 

Note there are certainly an infinite number of actions 
obeying the equivalence principle of Eq. ^3J Most of 
these will be non-local but there could well be more than 
one local action. The following hypothesis is made. If 
a local action obeying "strong" equivalence exists, then 
the measure p, r should fall exponentially, and the oper- 
ator D r should have exponentially falling coefficients in- 
side A r . D r is the best ultra-local approximation to the 
solution of the equivalence condition of Eq. ^| The coef- 
ficients of the action a long way from the boundary of the 
operator should also converge as the range is increased. 

The sequence of ultra-local actions D r was computed 
numerically by finding the minimum of ji r . The calcu- 
lations were performed for massless fermions. A short 
check demonstrated the localisation properties were bet- 
ter for massive fermions. 
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FIG. 1: The error function, fi r for lattice actions of range up 
to 106. 



A. Results 

A multi-dimensional Newton-Raphson solver was used, 
since both the slope and Hessian of p, r can be computed 
easily. The GNU multiple precision library (GMP) was 
used when numerical precision was required beyond 
64-bit native arithmetic. Some checks were made to test 
if the minimum in p, r was a global one. A range of dif- 
ferent starting values of the action parameters were used 
to seed the Newton-Raphson search and a simulated an- 
nealing algorithm was run to search for a minimum at 
short ranges. A number of local minima were found in 
many cases making it difficult to determine if the global 
minima was reached. This issue is discussed later. 



1. Two dimensions 

Fig. shows the dependence of p, r for the optimal 
action as a function of the finite range of the action, r. 
A clear signal for exponential fall-off is displayed: /i r , 
the discrepancy between D\D r and — □ falls by thirty 
decades as the action range is increased from b to 10&. 
The coefficients in q and p^ of the action Diq are pre- 
sented in Fig. [21 The on-lattice-axis and diagonal terms 
are presented. For the operator p^, on-axis refers to the 
terms in pi with off-set vector (a, 0) and those in j> 2 with 
off-set (0, a). Note that by symmetry, terms in p\ with 
off-set (0, a) vanish identically. An exponential fall-off 
over eighteen decades is observed, providing solid evi- 
dence for the existence of a local operator. At 64-bit 
machine precision, terms with ranges beyond about 66 
would have uncomputably small contributions to the ac- 
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FIG. 2: The coefficients in q and p M , the composite operators 
in the action Dio, as a function of the separation between the 
fields in the bilinear. The field separations are given using the 
usual 2-norm distance. 
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FIG. 4: The error function, fi r for four-dimensional lattice 
actions of range up to 5b. 
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FIG. 5: The coefficients in q and p M in the four- dimensional 
action D5 as a function of the separation between the fields 
in the bilinear 
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FIG. 3: The eigenvalues of the approximately equivalent 
Dirac operator, D10. The eigenvalues are scaled so the first 
doubler, with momentum (7r,0) has A = 1. 



tion of a plane wave state. Notice the terms in the two 
operators p^ and q are very similar in magnitude and 
have similar localisation range. Fig. |31 shows the eigen- 
value spectrum of the operator Din- The eigenvalues are 
purely real when all components of the operator p^ van- 
ish, it. at the doubling points. The real parts of the 
eigenvalues for the doublers are close to (for the prop- 
agating mode), 1 and v2- 



2. Four dimensions 

A similar construction was followed for the four- 
dimensional action. A sequence of operators on the four- 



dimensional lattice with spacing b = 2a was determined. 
The operators in the action again took the general form of 
Eqns. 1281 and 1241 with q transforming trivially under the 
four-dimensional rotation group and p^ transforming un- 
der the fundamental representation. The computational 
cost of the four-dimensional calculation means that only 
actions up to range r — 5 have been constructed. Fig. 
0] shows the fall-off of /i r in the sequence of actions. A 
rapid decay of p r as r is increased is again observed. No- 
tice also the fall-off accelerates at range r = 4, once the 
action increases beyond the bounds of the unit hypercube 
on the spacing-6 lattice. 

In Fig. the magnitude of the coefficients on the axis 
and diagonals of the operator D5 are displayed. A similar- 
pattern to the two-dimensional case is seen, with a five- 
decade decay over five lattice hops being observed. The 
data suggest a local operator exists in four dimensions as 
well as two. 
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IV. THE GINSPARG- WILSON RELATION 

The presence of the operator q in the definition of Eq. 
1191 means D does not anticommute with 75 , which would 
guarantee that the fermionic action is invariant under 
continuous chiral transformations. This is expected from 
the Nielsen-Ninomiya theorem. As is now well known, 
the theorem can be bypassed if one does not insist that 
the chiral transformations assume their canonical form 
on the lattice [Tol ITlj . In particular it was shown that 
the Ginsparg- Wilson relation, 



{ l5 ,D} = 2D l5 RD, 



(29) 



implies an exact symmetry of the fermionic action which 
may be regarded as a lattice form of an infinitesimal chi- 
ral rotation. 

For a Dirac operator obeying Eq. ^| three properties 
follow: 



D- 1 = 



75^75, 

-□ J, 

7m?V + g 
□ 



(30) 
(31) 

(32) 



The propagator D 1 then satisfies the following relation 



{^D- 1 } = 2i? 75 , 



with 



R - q 



(33) 



(34) 



The construction of Sec. IHII does not ensure the oper- 
ator R is local, and so there is no apparent lattice chiral 
symmetry. Eqn. 1341 does suggest the definition of an al- 
ternative sequence of actions, D^ GW ' which might lead 
to a Dirac operator obeying a (generalised) Ginsparg- 
Wilson relation. Consider an operator of range r with 
the form 



D (GW) 



7*«Pr 



DR. 



r— 1; 



where i? r _i is a local operator with finite range r 
Note this implies the scalar operator q r = DR. 



(35) 

- 1. 
1 has 



range r as before. If the limit of the sequence Dr W ' is 
a solution to Eqn. 1191 and p M and R remain local, then 
a local operator, equivalent to the staggered fermion and 
obeying a generalised Ginsparg- Wilson relation will be 
constructed. 

A simple consequence of Eq. [2!5]is then that chiral sym- 
metry is partly preserved, in particular the lagrangian 
L = ipDip is invariant under the local symmetry trans- 
formation: 



■><■ = --,( \-\rd), . 



5ip = ip[l- -DR I 75 



(36) 



(37) 
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FIG. 6: The error function, fi r for two-dimensional lattice 
operators constructed with the Ginsparg- Wilson constraint 
of Eqn. OH 
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FIG. 7: The coefficients in q and in the two-dimensional 
action D^ W ' as a function of the separation between the 
fields in the bilinear. 



results 



The operator sequence minimising /i r of Eq. 1271 was 
computed for the two-dimensional lattice. As before, the 
Newton solver was used for minimisation with the con- 
straint of Eqn. 1351 The problem of finding multiple local 
minima was observed and was more extreme than the 
initial construction of Sec. II I II The dependence of /i r 
for the operator obeying the Ginsparg- Wilson constraint 
is shown in Fig. Exponential decay of the discrep- 
ancy measurement, fi r as the range r is increased is seen 
again this time over twelve orders-of-magnitude. Fig. [7| 
shows the fall-off of the symmetry-breaking kernel for the 
optimised action D^f . An exponential fall-off of eight 
decades is observed. The terms in the operator q are thus 
exponentially localised. The derivative operator, also 
has local coefficients. The eigenvalues of the operator 
d[q W ^ are shown in Fig. along with those for the un- 
constrained two-dimensional action construction of Sec. 
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FIG. 8: The eigenvalues of the two operators, Dio (blue) and 
D[^ w ^ (red). The eigenvalues are scaled such that A = 1 for 
momentum (n,0). 



V. DISCUSSION 

A numerical construction of ultra-local, approximate 
actions can never prove the existence of a fermion with a 
well defined local action (unless the action is itself ultra- 
local) , but the calculations in this paper do present strong 
evidence for the existence of an equivalent, local the- 
ory to the one-flavour free staggered fermion. In two 
dimensions, the mis-match between the dispersion spec- 
trum of the ultra-local theory and the staggered fermion 
can be made as small as 10~ 30 with an action of range 
106. At this range, terms in the action are as small as 
10 -18 . The construction for four-dimensional theories is 
more difficult, but evidence for locality is seen here too. 
The numerical construction of the action was performed 
for massless fermions. Some short tests with massive 
fermions suggest the localisation properties of these ac- 
tions are better still; the massless fermion represents the 
hardest case to reproduce. 

The numerical search for a global minima of Eqn. 1271 
a non- linear function of the action parameters, lo p and 
ijj q is made difficult by the presence of local minima. It 
is difficult to find convincing evidence that the Newton- 
Raphson solver has found the global minimum for large 
actions, although searches using different starting points 
often converged to a common fixed point. An empiri- 
cal observation is important; the minima with smaller 
values of fi r tend to have better localisation properties 
(their coefficients fall faster). This implies if the searches 
have not found the global minima, these will represent 
better actions than those already uncovered, improving 
the construction rather than spoiling it. 



For the constrained construction to build an opera- 
tor obeying the (generalised) Ginsparg- Wilson relation, 
the localisation ranges were about twice that of the un- 
constrained construction and good evidence for exponen- 
tially local actions is seen. The problem of finding global 
minima seems to be exacerbated. Solutions to the stan- 
dard GW relation (with R = I) are now well known 
[l2L Il3j . It is important to recognise a solution with 
R = I is impossible for the equivalent theory; this can be 
seen easily by considering the doubler momenta, (7r, 0) 
and (7r,7r). For these two points, the eigenvalues of D 
must be 1 and \[2 respectively, while R — I would de- 
mand they were both unity. 

The evidence in this paper (and in Ref. J21.J) put 
staggered-fermion simulations on a more robust footing, 
but there remain many unanswered questions. For stag- 
gered fermion simulations to be correct descriptions of 
quantum field theory, one must demonstrate two postu- 
lates; firstly that a local path-integral representation of 
the fractional power of the staggered determinant exists 
and secondly the validity of calculations performed by 
assuming the propagator of this theory is related to the 
inverse of the full staggered matrix, Q . The work in 
this paper hints at the right question to address the first 
issue, but does not address the second point. In using 
Q _1 as the fermion propagator a clear problem arises. 
In four dimensions, too many pion operators can be con- 
structed for example. The residual symmetry of the stag- 
gered matrix ensures these states lie in mass-degenerate 
multiplets EHHEfl, but "taste-breaking" interactions 
split their masses. These splittings vanish in the con- 
tinuum limit. Recent work on the low-ener gy eigenvalue 
spectrum is beginning to resolve this issue flcl Il3. l20l| . 



Future work offers an optimistic possible outcome; if an 
effective operator can be constructed for QCD, it seems 
this operator might obey a Ginsparg- Wilson relation and 
the lattice physicist will have the best of both worlds: 
cheap dynamical-configuration generation algorithms us- 
ing the staggered formulation with a theoretically well 
defined action (possibly with an exact GW chiral sym- 
metry) to compute propagators. This would not be a 
"mixed action" simulation, where different valence and 
sea quark actions are employed. The first obstacle to 
extending the construction of Sec. 11111 to incorporate 
gauge interactions is that the operator Q*Q is not pro- 
portional to (I ® I) in Dirac-taste indices and so does 
not decompose directly into decoupled parts. This 
is the "strong" definition of equivalence required for the 
free theory, but a re-definition of \i r to measure viola- 
tions in "weak" equivalence can be made. This is under 
investigation and few conclusions about the success of 
this programme can be drawn. A number of difficult 
questions arise immediately, since the two theories would 
have different apparent symmetries. 
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VI. CONCLUSIONS 

In this paper, a local lattice Dirac operator whose de- 
terminant is identical to the free staggered-fermion de- 
terminant, and whose energy-momentum dispersion rela- 
tions are identical (although with different degeneracies) 
is described as the end-point of a sequence of actions of 
increasing, finite range. The first few ultra-local actions 
in the sequence are constructed numerically and conver- 
gence of the sequence is demonstrated in both two and 
four dimensions. 

The spectrum of the operator is free from doublers and 
its low-energy dynamics correctly describes the propaga- 
tion of free fermions up to corrections of C(a 2 p 2 ), a prop- 
erty it inherits from its staggered parent. The operator 
acts on a full Dirac spinor situated only on the sites of a 
blocked lattice with spacing b — 2a. 

A constraint is added to the construction to define an 



operator that obeys a generalised Ginsparg- Wilson rela- 
tion. In this action, the chiral symmetry breaking in the 
propagator is described by a local operator, diagonal in 
the spin index. This implies the existence of a fermion 
whose path integral is the same as that of one staggered 
flavour and with an exact chiral symmetry on the lattice. 
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